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ABSTRACT 

We present the SFI++ dataset, a homogeneously derived catalog of photometric and rotational 
properties and the Tully-Fisher distances and peculiar velocities derived from them. We make use 
of digital optical images, optical long-slit spectra, and global HI line profiles to extract parameters 
of relevance to disk scaling relations, incorporating several previously published datasets as well as 
a new photometric sample of some 2000 objects. According to the completeness of available redshift 
samples over the sky area, we exploit both a modified percolation algorithm and the Voronoi-Delaunay 
method to assign individual galaxies to groups as well as clusters, thereby reducing scatter introduced 
by local orbital motions. We also provide corrections to the peculiar velocities for both homogeneous 
and inhomogeneous Malmquist bias, making use of the 2MASS Redshift Survey density field to ap- 
proximate large scale structure. We summarize the sample selection criteria, corrections made to raw 
observational parameters, the grouping techniques, and our procedure for deriving peculiar velocities. 
The final SFI++ peculiar velocity catalog of 4861 field and cluster galaxies is large enough to permit 
the study not just of the global statistics of large scale flows but also of the details of the local velocity 
field. 

Subject headings: astronomical databases: miscellaneous — galaxies: distances and redshifts — galaxies: 
fundamental parameters — large-scale structure of the universe 



1. INTRODUCTION 

The recessional velocities of galaxies exhibit deviations 
from smooth Hubble flow which are induced by inho- 
mogeneities in the matter distribution of the universe. 
These deviations from smooth Hubble flow are known as 
peculiar velocities. Because the peculiar velocity field is 
determined by the scale and amplitude of the inhomo- 
geneities, the measurement of galaxy peculiar velocities 
can be used to trace the distribution of the mass, both 
dark and light. 

The galaxy peculiar velocity field is thus a powerful 
tool. It allows us to map the large scale structures of 
the local universe, independent of the distribution of lu- 
minous matter. And because the large scale structure 
depends on cosmological parameters, it can also be used 
to derive cosmological parameters, including the cosmo- 
logical matter density, rate of growth of structure, and 
Hubble expansion rate. 

The measurement of peculiar velocities is unavoidably 
intertwined with the measurement of distances, in that 
redshift-independent distances are needed in combina- 
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tion with redshifts to extract peculiar velocities. Given 
a galaxy's observed recessional velocity cz and distance 
r, as measured by a redshift-independent distance in- 
dicator, one can infer its line-of-sight peculiar velocity 
according to 

Vpec =cz- H Q r (1) 

Here and in the remainder of this work, in order to avoid 
confusion between recessional velocities and peculiar ve- 
locities, recessional velocities are designated by cz, and 
peculiar velocities by v. Both quantities will be accom- 
panied by a subscript where appropriate. All redshifts 
are assumed to be measured in the cosmic microwave 
background frame unless otherwise noted. 

Most of the largest peculiar velocity surveys have em- 
ployed secondary distance indicators — in particular the 
Fundamental Plane (Djorgovski & Davis 1987, Dressier 
et al. 1987; hereafter "FP") and Tully-Fisher relations 
(Tully & Fisher 1977; hereafter, "TF relation"). The 
former expresses the luminosity of an elliptical galaxy as 
a power law function of its radius and velocity disper- 
sion, while the latter expresses the luminosity of a spiral 
galaxy as a power law function of its rotational velocity. 
The earliest surveys included peculiar velocity measure- 
ments for of order ~ 1000 galaxies at most. Many of 
these early surveys were concatenated together into the 
Mark III catalog (Willick et al. 1995, 1996, 1997). The 
earliest catalogs containing more than 1000 galaxies in- 
clude the Spiral Field I-band (SFI; Giovanelli et al. 1994, 
1995, Haynes et al. 1999a, 1999b), Spiral Cluster I-band 
(SCI; Giovanelli et al. 1997a, 1997b), Spiral Cluster I- 
band 2 (SC2; Dale et al. 1999a, 1999b), ENEAR (da 
Costa et al. 2000, Bernardi et al. 2002), and EFAR (Col- 
less et al. 2000). While these surveys have been useful 
for illuminating the global features of large scale motions 
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in the universe, there are still some inconsistent results 
relating to fundamental aspects of the motions, such as 
the scale of the largest flows and the value of (3 (see the 
review papers Dekel 1999 and Zaroubi 2001). Larger 
samples are needed to resolve these discrepancies, and to 
investigate the details of the large scale flows. 

We present here the SFI++, one of the first of a new 
generation of peculiar velocity surveys of ~ 5000 galax- 
ies or more. The catalog includes I-band TF parame- 
ters for 4861 field and cluster galaxies. The only pecu- 
liar velocity catalog presently in existence to exceed the 
size of SFlH — \~ is the Kinematics of the Local Universe 
catalog (KLUN; Theureau 1998, Theureau et al. 1998, 
and references therein), which consists of 6600 galaxies 
with apparent B magnitudes and HI redshifts and line 
widths. The B-band TF relation has significantly more 
scatter than the I-band relation, resulting in larger dis- 
tance errors. This is because both Galactic and internal 
extinction are more significant at B-band, and because 
the stars that supply the light are largely confined to 
star-forming regions — not distributed as smoothly as the 
stars in I-band — making accurate disk inclinations more 
difficult to measure. Additionally, a significant fraction 
of the data in KLUN are from the literature, and are 
at somewhat low velocity resolution. In contrast, the 
SFI++ datasets contain both new photometry derived 
from isophotal fitting of images and digital optical rota- 
tion curves (ORCs) and HI profiles for which methods of 
measuring widths have been derived specifically for the 
purpose of recovering the circular rotational velocity. 

Observations are currently underway on the KLUN+ 
TF observing program, which aims to expand KLUN 
to include 20,000 spiral field galaxies with Nancay ra- 
dio telescope HI spectra, and added B and I band pho- 
tometry (Theureau et al. 2005). Other large peculiar ve- 
locity surveys with observations in progress include the 
NOAO Fundamental Plane Survey (Smith et al. 2004), 
which aims to provide FP measurements for 4000 early- 
type galaxies in 100 x-ray selected clusters within 200 
Mpc/h, the 2MASS Peculiar Velocity Survey (Masters et 
al. 2005), which is planned to include TF measurements 
for more than 5000 of the brightest inclined spirals in 
the 2MASS Redshift Survey (2MRS; Huchra et al. 2005) 
and could be extended to include many more such ob- 
jects, and the 6dF Galaxy Survey (Jones et al. 2004), 
which is planned to include optical photometry and line 
widths, providing FP distances for 15,000 galaxies. 

The TF relation requires three observational compo- 
nents for determining distances (and thus peculiar ve- 
locities): galaxy systemic velocities (redshifts), appar- 
ent magnitudes, and rotational velocity widths (which 
also require disk inclination estimates) . Our group main- 
tains large databases of 21 cm spectral line parameters 
(Springob et al. 2005, hereafter S05), optical rotation 
curve parameters (Catinella, Haynes, & Giovanelli 2005) 
and I-band photometric parameters based on observa- 
tions made over the course of the last 20+ years. We 
have used subsets of these data for several past TF stud- 
ies. However, as explained in the aforementioned pa- 
pers, much of this observational dataset has been newly 
reprocessed. In addition, a large number of photomet- 
ric observational parameters have never been published 
until now. We have subsequently synthesized the en- 
tirety of these data into a new TF catalog, the SFI++, 



which we present here. Because the new radio and op- 
tical spectroscopy data have already been presented by 
S05 and Catinella, Haynes, & Giovanelli (2005) respec- 
tively, we have not replicated the exhaustive compilation 
of spectroscopic parameters here. Instead, we include 
a handful of parameters most directly applicable to the 
computation of TF distances. In contrast, there is no 
corresponding publication of our newly reprocessed pho- 
tometric data, though an earlier generation of the catalog 
was presented by Haynes et al. (1999a, hereafter H99), so 
we provide here additional photometric parameters that 
are not directly used in the calculation of TF distances. 

If one is to use the TF relation to compute peculiar 
velocities, one must first calibrate the power law with a 
template relation. The SCI catalog was the first attempt 
by our group to derive such a template relation. The 
derivation of the SCI template was presented by Gio- 
vanelli et al. (1997b), while the SCI dataset itself was 
presented by Giovanelli et al. (1997a, hereafter G97). 
We have similarly divided our presentation of SFI++, 
such that the derivation of our template relation is pre- 
sented in Masters et al. (2006; hereafter Paper I), while 
this work is devoted to the discussion and presentation 
of the catalog itself. As explained in Paper I, a new 
template relation is required here because of the signif- 
icantly larger data sample and the revised corrections 
to raw data. But unlike G97, which just presented the 
cluster data, we present the data for both clusters and 
field galaxies in this work. In future papers, including 
Springob et al. (2007, in prep.) and Masters et al. (2007, 
in prep.), we will use this dataset to investigate the local 
peculiar velocity field. 

In Section 2, we review the observational selection cri- 
teria that was used for each of the individual observing 
campaigns that define the SFI++ sample. In Section 3, 
we discuss the process by which raw observational pa- 
rameters are corrected to produce physically meaningful 
estimates of the galaxies' luminosities and rotational ve- 
locities. In Section 4, we present additional photometric 
parameters for the galaxies, separate from those directly 
used in the derivation of peculiar velocities. In Section 5, 
we describe the cluster template we have derived in Paper 
I, which we use to calibrate the TF relation. In Section 
6, we describe the origin of the group assignments that 
we have adopted. In Section 7, we explain the derivation 
of peculiar velocities for individual galaxies and galaxy 
groups. The SFI++ catalog itself is presented in Section 
8, and catalog properties are discussed in Section 9. A 
brief summary concludes in Section 10. 

2. SELECTION OF SAMPLE 

The complete photometric and spectroscopic data sam- 
ple is composed of the previously published SFI, SCI, and 
SC2 datasets, and the until now unpublished Spiral Field 
I-band 2 (SF2) sample. In each of these cases, prospec- 
tive targets were selected from among spirals included 
in our private database, referred to as the Arecibo Gen- 
eral Catalog (AGC). We also include the TF samples of 
Mathewson, Ford, & Buchhorn (1992) and Mathcwson & 
Ford (1996), but reprocessed to extract parameters using 
our methodology. 

As explained in these references: the HI spectroscopy 
observations were made with the 305 m Arecibo telescope 
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of the National Astronomy and Ionosphere Center 7 , the 
late 91 m and 42 m Green Bank telescopes of the Na- 
tional Radio Astronomy Observatory 8 , the Nancay ra- 
dio telescope of the Observatory of Paris, and the Ef- 
felsberg 100 m telescope of the Max Planck Institut fiir 
Radioastronomie; the optical spectroscopic observations 
were made with the 2.3 m telescope at Siding Spring 
Observatory, the Hale 5 m telescope at Palomar Obser- 
vatory 9 , and the Cerro Tololo Inter- American Observa- 
tory (CTIO) 4 m telescope; and the optical photometry 
was done with the 1 m and 3.9 m telescopes at Siding 
Spring Observatory, the Kitt Peak National Observatory 
(KPNO) and CTIO 10 0.9 m telescopes, and the 1.3 m 
McGraw-Hill telescope of the Michigan-Dartmouth-MIT 
(MDM) Observatory 11 . The details of the selection, ob- 
servations, and data reduction are covered in each of the 
aforementioned papers. However, below we recap some 
of the main features of the sample selection for each of 
these individual projects. 

The SCI was a compilation of TF measurements for 24 
clusters, compiled for the determination of a template 
TF relation and the determination of the motions of the 
clusters themselves. The selection of the clusters used in 
SCI is described in G97 Section 2.3. The clusters all have 
mean velocities less than 10,000 km s _1 when measured 
in the CMB frame. They were chosen so as to span a 
large range in richness, and provide as much balance as 
possible among different parts of the sky. As explained 
in G97, the SCI itself includes some data from the lit- 
erature, including Pierce & Tully (1988), Han (1992), 
and Han & Mould (1992). However, the raw observa- 
tional parameters from these supplementary data were 
corrected using the same algorithms that were applied 
to the data from our own group. All of the data were 
processed in the same way. 

The SC2 was a TF compilation of clusters, conducted 
to improve the quality of the I band TF relation and 
to determine the redshift depth of the structures most 
responsible for the reflex motion of the Local Group with 
respect to the CMB. 5-15 TF measurements per cluster 
were obtained for an all-sky sample of 52 clusters with 
recessional velocities 5000 < cz < 25, 000 km s _1 . 

The SFI was comprised of a TF sample of 2000 field 
galaxies limited to cz < 7500 km s _1 (Local Group 
frame), blue magnitude tub < 14.5, and line width > 100 
km s _1 . The SFI also had redshift dependent upper and 
lower optical diameter limits, to minimize the variations 
in the number of objects observed per unit redshift that 
show up in all flux or diameter limited catalogs. The di- 
ameter limits were 2.5' < a < 5.0' for cz < 3000 km s _1 , 

7 The National Astronomy and Ionosphere Center is operated 
by Cornell University under a management agreement with the 
National Science Foundation. 
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the National Science Foundation. 

9 The Hale Telescope is operated by the California Institute of 
Technology under a cooperative agreement with Cornell University 
and the Jet Propulsion Laboratory. 

10 KPNO and CTIO are operated by Associated Universities for 
Research in Astronomy, under a cooperative agreement with the 
National Science Foundation. 

11 The MDM Observatory was jointly operated by the University 
of Michigan, Dartmouth College, and the Massachusetts Institute 
of Technology on Kitt Peak mountain, Arizona. 



1.5' < a < 5.0' for 3000 < cz < 5000 km s" 1 , and 
1.3' < a < 5.0' for 5000 < cz < 7500 km s" 1 . 

The SF2 program was intended to obtain photometry 
for objects either with existing HI or optical spectroscopy 
in our existing database at the time of the observations 
and to target the region —15 < decl. < +35° to a depth 
of cz < 10,000 km s" 1 and optical diameter a > 0.9'. 
I band images were obtained at the KPNO 0.9 m tele- 
scope. Approximately 2300 23' x 23' fields were observed, 
of which roughly 1900 yielded high quality photometry. 
Roughly 500 of these fields were used as part of SC2, with 
the remaining fields used for SF2. The targeted fields 
contained at least one good TF candidate (undisturbed, 
inclined spiral with cz < 10000 km s _1 ) for which a 
rotational width was already available to us from either 
HI line or optical long slit spectroscopy at the time the 
photometry was obtained. Moreover, since additional 
spectroscopic observations were possible, the 23' x 23' 
fields were centered to maximize the number of poten- 
tial additional TF candidates, regardless of the status of 
rotational width measurements. As a result, the number 
of galaxies for which photometry is available exceeds the 
number with rotational velocities. Additional factors re- 
lated to the allocation of telescope time, the eventual 
weather conditions and the practicalities of observing 
also had an impact on the final sample. For example, 
more observing time was allocated in the fall than in the 
spring, and the fall time had better weather. Also, while 
the original plan was to observe only fields with \b\ > 20°, 
some fields closer to the Galactic plane were included to 
fill parts of nights when higher latitudes where not ac- 
cessible. 

The Mathcwson, Ford, & Buchhorn (1992) dataset was 
a compilation of TF measurements for 1355 Sb-Sd galax- 
ies in the southern hemisphere, diameters a > 1.7', in- 
clinations > 40°, and Galactic latitude \b\ > 11°. Most 
of the objects had systemic velocities of less than 7000 
km s _1 . However, in the "Great Attractor" region, some 
higher redshift galaxies, drawn from the redshift survey 
of Dressier (1988), were included. The Mathewson & 
Ford (1996) dataset was a compilation of TF measure- 
ments for an additional 920 Sb-Sc galaxies selected from 
the ESO-Uppsala Survey of the ESO(B) Atlas (Lauberts 
1982, hereafter ESO), with diameters 1.0' < a < 1.6', 
systemic velocities between 4000 and 14,000 km s _1 , 
and the same inclination and Galactic latitude limits as 
the Mathewson, Ford, & Buchhorn (1992) sample. An 
additional 172 Uppsala General Catalog (Nilson 1973, 
hereafter UGC) galaxies were observed in the region 
250 < Z < 360°, 45 < b < 80°. Both the photomet- 
ric and optical spectroscopic datasets from Mathewson, 
Ford, & Buchhorn (1992) and Mathewson & Ford (1996) 
were made available to us and have been reprocessed us- 
ing our own algorithms to achieve greater homogeneity 
(Haynes et al. 1999a). 

The TF sample presented here, which we refer to as 
SFI++, is the union of each of these datasets. We in- 
clude only those observations deemed to be of 'high qual- 
ity', so that, for example, HI spectral profiles assigned 
the quality index 'G' as defined by S05 Section 4 are in- 
cluded, but all other HI spectra are not. Some galaxies 
have photometry, but not spectroscopy (or at least not 
of high quality) . Others have spectroscopy but not pho- 
tometry. Such objects are excluded from consideration 
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in the peculiar velocity catalog. However, we do include 
galaxies without high quality spectroscopy in a separate 
photometric compilation that will be presented in Sec- 
tion 4. 

There are also some objects with multiple photometric 
or spectroscopic observations. In such cases, we choose 
the observations deemed to be of highest quality, and 
discard the others. Preference is given to HI spectroscopy 
over optical spectroscopy, and to our own observations 
over those from the supplementary datasets. 

The final sample includes 4861 galaxies for which we 
have good TF data, of which 807 are members of the 
template clusters, which we discuss in Section 5. 

3. TULLY-FISHER PARAMETERS 

3.1. Photometry 

3.1.1. I band fluxes 

I-band optical photometric images were reduced as 
described by H99. Corrected I-band apparent magni- 
tudes and errors are then computed as described in that 
work. However, while our previously published photom- 
etry used Galactic extinction corrections derived from 
Burstein & Heiles (1978), we now use updated values 
taken from the Diffuse Infrared Background Experiment 
on the Cosmic Background Explorer Satellite (Schlegel, 
Finkbcincr, & Davis 1998). For estimating extinction 
from the target galaxy itself, we use Equation 27 from 
Giovanelli et al. (1994): AM = -jlog(a/b), where AM 
is the extinction in magnitudes, a and b are the observed 
semimajor and semiminor axes of the galaxy obtained 
from ispohote fitting, and 7 is a quantity that depends 
on the galaxy's inferred absolute magnitude, as described 
by that paper. 7's dependence on absolute magnitude 
Mi is shown in Giovanelli et al. (1995) Figure 7c. The 
exact functional form, which we adopt here, is 



7 = 0.5 for Mi > -19.1 



(2) 



7 = l-0.417(M/ + 20.3) for -20.3 < Mj < -19.1 (3) 
7 = 1.0 for - 21.8 < Mi < -20.3 (4) 



7 = 1.35 - 0.35(M/ + 22.8) for - 22.7 < M/ < -21. 



7 = 1.30 for Mi < -22.7 



(5) 
(6) 



We note that these relations were derived for a sample 
that mainly consisted of Sbc and Sc galaxies, whereas the 
sample presented in this paper includes a greater diver- 
sity of spiral subclasses, most notably some earlier mor- 
phological types. One would expect that, since we use ax- 
ial ratio as a tracer of inclination, and earlier morpholog- 
ical types tend to have thicker disks than later types, we 
may be underestimating the extinction in earlier types. 
However, this effect should be at least partially offset by 
the fact that later types are likely to have more dust. In 
any case, because we have derived separate morpholog- 
ical corrections to the TF template relation in Paper I, 
any underestimate of the extent to which 7 varies with 
morphological type should manifest itself as an increase 
in the TF scatter for the earlier types. There should 
be no systematic effect, as the galaxies of any particular 
inclination are randomly distributed across the sky. 



3.1.2. Inclinations 

Once the photometric images are reduced, inclinations 
are then computed from the data as in G97 Section 4: 
That is, from the measured semimajor and semiminor 
axes a and 6, we obtain the ellipticity e = 1 — 6/aasa 
function of the distance r = a/2 from the center of the 
galaxy. A range of radial distances is then chosen over 
which the disk appears to be exponential, and we then 
obtain a mean value of the ellipticity e. This elliptic- 
ity is then corrected for the smearing effects of seeing, 
according to G97 Equation 2: 



1 



'(1-e) 2 - V 2 



1 - VJ 2 



(7) 



where V is 2.5 times the size of the seeing disk for 
the observation, divided by the major axis at the 23.5 
mag / arcsec 2 I-band isophote. The inclination angle i is 
then given by G97 Equation 3: 



(cosi) 2 



(1 e corr ) 



1 



9o 2 



(8) 



where qo is the intrinsic axial ratio of the disk, taken to 
be 0.13 for galaxies of morphological type Sbc and later, 
and 0.20 for galaxies of earlier types. 

3.2. Spectroscopy 
3.2.1. 21 cm velocity widths 

Raw 21 cm line velocity widths were corrected for in- 
strumental and noise effects as described in S05, Sec- 
tion 3.2.2. That is, the widths are measured by taking 
the difference of velocities between the two midpoints 
of polynomials fit to either side of the line profile. The 
midpoints arc defined as the point at which the flux (as 
given by the polynomial) is 50% of the peak minus rms 
value. In the vast majority of cases, we used a first order 
polynomial. But in rare cases, a second order polyno- 
mial was used. Instrumental and noise corrections are 
applied, as is a redshift correction, all as in S05, Section 
3.2.2. 6.5 km s" 1 was then subtracted from this width as 
described by S05, Section 3.2.3, to produce a turbulence 
corrected width, which S05 refers to as W c ,t- This value 
was then corrected for line-of-sight projection effects by 
dividing by the sine of the inclination angle i, where the 
inclination angle was derived as in Section 3.1.2. 

In summary then, the corrected HI widths W21 are 
given by 



21 



W obs , 2 i - A s 
1 + z 



1 

sim 



(9) 



where W b s ,2i is the observed width, and A s and A t 
are the instrumental and turbulence corrections respec- 
tively. (This expression only differs from G97 Equation 
5 in that we subtract the turbulence correction linearly 
rather than quadratically.) As in S05, we use A t = 6.5 
km s _1 , and A s = 2AvX, where Av is the spectrometer 
channel separation in km s" 1 and A is a function of the 
SNR and type of smoothing as described in S05 Section 
3.2.2. Sec S05 Table 2 for the precise dependence of A 
on SNR and smoothing type. 

Errors on W21 are computed exactly as in G97 Equa- 
tion 7, with all terms defined as in that paper. That 
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is, the errors are computed as the sum in quadrature 
of the errors on the observed width, instrumental and 
noise corrections, turbulence correction, and inclination 
correction. 

3.2.2. Optical velocity widths 

All observed optical velocity widths have been ex- 
tracted by fitting a function to the folded Ha ORCs. 
We use the parametric Polyex model first described by 
Giovanelli & Haynes (2002). As described by Catinclla, 
Haynes, & Giovanelli (2005), this model has the func- 
tional form for the circular rotational velocity Vpg at a 
distance r from the galaxy's center 

V PE (r) = V (l - e- r ^ E )(l + ar/r PE ) (10) 

Vo and tpe are, respectively, the circular velocity am- 
plitude and exponential scale of the inner region of the 
galaxy, a is the slope of the outer part of the ORC. Ob- 
served widths are measured using the value of this func- 
tion at r opt , an optical radius containing 83% of the total 
light of the galaxy. Corrected widths are then computed 
from the observed widths, as described by Catinclla, 
Haynes, & Giovanelli (2005), using the same cosmolog- 
ical and inclination corrections that are used for 21 cm 
line widths. 

Catinella, Haynes, & Giovanelli (2007) shows that 
there are slight systematic differences between widths 
measured from HI spectroscopy and widths measured 
from optical rotation curves. These differences depend 
on the relative extent of the Ha (r max ) as compared to 
the total optical extent of the galaxy, and on the slope 
of the ORC at the optical radius. We thus correct each 
of the ORC widths using the relation: 

W21/W0RC = 0.899 + 0.188r max /r opt for rising ORCs 

(11) 

W 21 /W rc = 1-075 - 0M3r max /r opt for flat ORCs 

(12) 

as derived in Catinella, Haynes, & Giovanelli (2007). 
Flat ORCs are defined as those for which the gra- 
dient of the rotation curve at r opt is less than 0.5 
km s _1 arcscc -1 , while rising ORCs have gradients 
greater than 0.5 km s _1 arcsec -1 . 

Both optical and radio widths for which all corrections, 
including inclination angle corrections, have been applied 
are hereafter designated by 'Wtf'- 

3.2.3. Recessional velocities 

For each galaxy, the recessional velocity cz is taken to 
be the midpoint of the spectral line profile, either HI pro- 
file or ORC, as explained in S05 and Catinella, Haynes, 
& Giovanelli (2005) respectively. As in our previous TF 
studies, we neglect redshift errors in computing TF dis- 
tances, as they are typically less than 1%. 

4. ADDITIONAL PHOTOMETRIC DATA 

We provide additional photometric parameters, of 
varying levels of applicability to the computation of TF 
distances, for an overlapping sample of 5254 galaxies. 
Excluded from this overlapping photometric sample are 
the SC2 galaxies and much of the data we have taken 



from the literature, for which we are missing these ad- 
ditional photometric parameters. However, we include 
galaxies for which we do not have high quality width 
measurements, which are thus not in the TF sample. 

The scale length, surface magnitude, ellipticity, posi- 
tion angle, and observed magnitude parameters are all 
defined exactly as in H99, Sections 2 and 3 (with all 
of the same corrections applied), which itself builds on 
the data reduction routines developed in Giovanelli et 
al. (1994). Note that our definition of the ellipticity is 
e = 1 — b/a, where a and b are the major and minor 
axes, respectively. (In this work, unlike in H99, we use e 
to denote ellipticity, as e is used to denote errors.) 

We provide these photometric parameters in Table 1. 
All galaxies for which we have these additional param- 
eters are included here. The format of the table is as 
follows: 

Column (1). — Entry number in the UGC, where appli- 
cable, or else in the AGC. 

Column (2). — NGC or IC designation, or other name, 
typically from the Catalog of Galaxies and Clusters of 
Galaxies (Zwicky et al. 1961-1968), ESO, or the Mor- 
phological Catalog of Galaxies (Vorontsov-Velyaminov & 
Arhipova 1968). Where used, the designation in the lat- 
ter is abbreviated to eight characters. 

Columns (3) and (4). — Right ascension (in hh mm 
ss.s format) and declination (in dd mm ss format) in 
J2000.0 epoch either from the NASA/IPAC Extragalac- 
tic Database (NED) 12 or measured by us on the POSS- 
I 13 . Typically, the listed positions have < 5"accuracy. 

Column (5). — The morphological type code following 
the RC3 system. Classification comes either from the 
UGC or ESO catalogs, or from our own visual examina- 
tion of the POSS-I prints. 

Column (6). — The isophotal radius measured at an I- 
band surface magnitude of 23.5 mag/arcsec 2 , ^3.5, in 
arcseconds. 

Column (7). — The optical radius, r^3L, in arcseconds, 
derived from the image and corresponding to the radius 
encompassing 83% of the light. 

Column (8). — The surface magnitude at the outer- 
most detected isophote in the I-band image, fi ut, in 
mag/arcsec 2 . 

Column (9). — The disk ellipticity, e, corrected for see- 
ing following Giovanelli et al. (1994). 

Column (10). — The estimated error on the disk ellip- 
ticity, e e . 

Column (11). — The mean disk position angle, PA, in 
degrees, computed by averaging the position angles from 
a series of isophotal fits of the galaxy. 

Column (12). — The estimated error on the position an- 
gle, e PA - 

Column (13). — The observed I-band magnitude, m b s , 
extrapolated to 8 disk scale lengths, before the extinction 
and face-on corrections have been applied that produce 
the values used for TF. 

Column (14)- — The measurement error on the appar- 

12 The NASA/IPAC Extragalactic Database is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, un- 
der contract with the National Aeronautics and Space Administra- 
tion. 

13 The National Geographic Society - Palomar Observatory Sky 
Atlas (POSS-I) was made by the California Institute of Technology 
with grants from the National Geographic Society. 
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Fig. 1. — Distribution of (a) cllipticity, in bins of width 0.05, 
(b) surface magnitude at the outermost isophotc, in bins of width 
0.1 magnitudes per square arcsecond, (c) isophotal radius at which 
the surface magnitude is 23.5 mag/arcscc 2 , in bins of width 5 arc- 
seconds, (d) observed I-band magnitude, in bins of width 0.2 mag- 
nitudes, (e) optical radius encompassing 83% of the I-band light, 
in bins of width 5 arcscconds, and (f) measurement error on the 
apparent I-band magnitude, in bins of width 0.01 magnitudes, for 
the 5254 galaxies for which we have these parameters. Note the 
difference in vertical scale between panels e and /, in contrast to 
the identical vertical scales for the other adjacent pairs. Figure 2 
of H99 shows the distribution of the same parameters, albeit with 
some slightly different nomenclature, for the Sc sample that existed 
at that time. 

cnt I-band magnitude, e m . 

Column (15). — A code indicating whether the galaxy is 
included in the TF sample. Those objects marked with 
a '*' are not included in the TF sample, as we do not 
have high quality width measurements for those galaxies. 
Any object not marked with a '*' should appear in cither 
Table 2 or Table 4. 

Figure 1 shows the distribution of each of the param- 
eters given in Table 1: e , r 23 .5, r 83L , [i ouU m obs , and 
e m , for the galaxies for which we have these additional 
parameters. H99 shows the histograms for exactly the 
same parameters for the corresponding dataset presented 
in that paper. Like the H99 sample, this sample is not 
magnitude limited. In fact, the distributions of r^.bi 
r%3L, and \i ou t are all remarkably similarly to the H99 
distributions. There do appear to be subtle differences in 
the ellipticity and magnitude error distributions however. 
The typical apparent magnitude measurement error is 
~ 0.03, as opposed to ~ 0.04 in H99. Also, the elliptic- 
ities of both samples peaks at about 0.75-0.80, but the 
dropoff towards low inclination galaxies is more gradual 
in the sample presented here. This more gradual dropoff 
is probably due to the inclusion of earlier morphological 
types, which have a larger intrinsic axial ratio, and to 
the fact that, unlike the SFI, SFI++ has no strict axial 
ratio selection criteria. 

5. SFI++ CLUSTER TEMPLATE 
We can express the TF relation as 



L = W° 



(13) 



where L is the luminosity of the galaxy, V max is the max- 



imum rotational velocity of the galaxy, and fc and a are 
constants that need to be calibrated with a TF template 
relation. As galaxy clusters provide a large number of ob- 
jects that are located at a common distance, they are well 
suited for the construction of such a template. The inde- 
pendent distance estimates of several galaxies in a clus- 
ter also provides a much more accurate determination of 
the cluster distance than is available for a single galaxy, 
providing us with "hard points" in the peculiar velocity 
field. While some early TF studies depended on tem- 
plates derived from a single cluster (Pierce & Tully 1988; 
Mathewson, Ford, & Buchhorn 1992), this approach of- 
fers significant drawbacks, as outlined by G97. Giovanelli 
et al. (1997b) developed the "basket of clusters" method 
to alleviate the bulk of these problems. In Paper I we 
have applied this technique to the larger SFI++ template 
sample to provide a new I-band TF template specifically 
designed for use with SFIH — \- galaxies. 

There are 86 different clusters for which we have data 
in either SCI or SC2. The cluster selection and galaxy 
cluster assignments are described by those papers. We 
will just mention briefly that the vast majority of these 
are Abell clusters, though there are a few exceptions 
in the SCI. As described by the aforementioned au- 
thors, cluster assignments were made by eyeball esti- 
mates based on the spatial and redshift distribution of 
galaxies near the cluster centers. 

All of these clusters which include at least seven mem- 
bers in SFI++ and average systemic velocities less than 
10,000 km s _1 in the CMB frame were considered for 
inclusion in the template. While they were already used 
as template clusters for either SCI or SC2, the much 
larger data sample that we have with SFI++ compelled 
us to reexamine the assignments of each of the galax- 
ies in the vicinity of these prospective template clusters. 
The redshift space positions of each of the galaxies in all 
prospective clusters was examined by eye, just as in G97 
and Dale et al. (1999a, 1999b). Following the scheme 
of G97, we differentiate between objects believed to be- 
long to the template cluster (the in sample), and objects 
with velocities very close to the systemic velocity of the 
cluster but spatially removed from their center so that a 
firm membership assignment cannot be made (the in+ 
sample). As much as possible, we have used the same 
thresholds for determining membership that were used 
for SCI and SC2. 

After each prospective template cluster was examined, 
we discarded a few clusters because they either had too 
few in objects with the new cluster assignments, or they 
were close enough to the galactic plane that we had con- 
cerns about differential Galactic extinction within the 
cluster. We have also omitted the Virgo cluster from the 
template, as uncertain assignments to its various sub- 
groups can create significant complications. The remain- 
ing 31 clusters that match our criteria, of which 22 are 
from SCI and 9 are from SC2, comprise the template. 
The resulting TF template relation, as derived in Pa- 
per I, expresses the corrected I band magnitude Mj as a 
function of the corrected rotational width Wtf as 

Mi = - 7.85 [log( W TF ) - 2.5] - 20.85 + 51og(/i). (14) 

Using the terminology of Paper I, this template uses 
the bivariate fit to the in+ sample, with the width- 
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dependent morphological correction applied. That is, 
because of the fact that the TF relation has some depen- 
dence on morphological type, Paper I uses the template 
relation of Equation 14 for galaxies with morphological 
types of Sbc and later only, while the magnitudes of ear- 
lier type galaxies are corrected by: 

• SO/Sa/Sab: -0.32 - 0.9[log(HW) - 2.5] mag 

• Sb: -0.10 - 0.9[log(W T F) - 2.5] mag 

These corrections were derived in Paper I by fitting the 
offset and difference in slope between the TF relation of 
the later type galaxies, and the TF relations of the earlier 
type galaxies. As Paper I explains, if the shallower slope 
of the TF relation for earlier type galaxies is intrinsic, 
it suggests that they have less concentrated halos than 
later types at a given rotational velocity. However, the 
difference in slope could also be due to differing levels of 
incompleteness for earlier and later type spirals. 

As with the SCI and SC2, the template fitting proce- 
dure also produces estimates of the peculiar velocities of 
the template clusters themselves. This is also described 
in Paper I. 

G97 explains the rationale behind the selection of the 
SCI/SC2 clusters, and the same principles apply here as 
well. It is desirable to have 'spatial balance', such that 
the clusters are distributed evenly across the sky, and 
across a range of redshifts. This increases the probability 
that the 'basket of clusters' is at rest with respect to 
the CMB. Including clusters across a range of redshifts 
also allows us to sample a large dynamic range of TF 
parameters. 

6. NONTEMPLATE GROUPS 

We also exploit the group assignments of non-template 
galaxies so that the measured distances of each galaxy 
in a group can be averaged, and the TF scatter on these 
points can be reduced by a factor ~ >/N for a group with 
N galaxies with TF measurements. Group assignments 
for the non-template galaxies are discussed below. 

Different authors have used a wide variety of group 
identification methods, the most common being the hi- 
erarchical (e.g., Tully 1987) and percolation (also known 
as "friends-of- friends" ; e.g., Huchra & Geller 1982) meth- 
ods. The fact that both the AGC and the HI archive are 
inhomogeneous catalogs with no easily describable selec- 
tion criteria creates significant complications for the as- 
signment of group identifications. However, we do have 
one significant advantage in this endeavor, in that we are 
not interested in studying the groups as groups. We only 
need the groups to provide us with more accurate dis- 
tance measurements. We are thus justified in using a het- 
erogeneous mixture of group identifications and galaxy 
group assignments, which we describe below. 

There are three different sources from which we have 
drawn group assignments. The first two, the SCI/SC2 
clusters (discussed in Section 5) and Nearby Optical 
Galaxy (NOG) P2 groups and their extensions, have been 
combined into a single group catalog, which was used in 
the determination of the HIMF by Springob, Haynes, & 
Giovanclli (2005). The third source of group identifica- 
tions, which we will discuss in Section 6.2, are Voronoi- 
Delaunay method (VDM) groups (Marinoni et al. 2002). 

Note: Particularly dense concentrations of galaxies are 
commonly referred to as clusters. Dense clusters have 
markedly different properties from loose groups of galax- 



ies, and the terms 'group' and 'cluster' are commonly 
defined in such a way that the two classifications are con- 
sidered to be mutually exclusive. However, group iden- 
tification algorithms are usually designed in such a way 
as to identify clusters as well. In the remainder of this 
work, except where stated otherwise, we define 'group' 
broadly, so as to include clusters as well. However, since 
the template clusters are all actual clusters, we continue 
to refer to them as 'clusters'. 

6.1. NOG groups 

The basic idea involved in a percolation group-finding 
algorithm is presented by Huchra & Geller (1982). In 
brief: for each galaxy in the catalog in question, one 
searches for neighboring galaxies within a given search 
radius in rcdshift space in both the transverse and radial 
directions. Any neighboring galaxies found to be within 
Dl in the transverse direction and Vl in redshift is as- 
signed to the same group as the original galaxy. Any 
such neighbors of the neighboring galaxies are also as- 
signed to the same group, and so on. Dl and Vl are 
referred to as 'linking parameters', and their values may 
be fine-tuned to identify systems of a particular number 
density contrast. 

In any flux or diameter limited survey, the density of 
detected objects drops off with increasing redshift. Thus, 
in order for a percolation algorithm to identify groups of 
the same density contrast at different redshifts, the link- 
ing parameters are usually allowed to vary with redshift. 

NOG is a complete, distance-limited (cz < 6000 
km s _1 , Local Group frame) and magnitude limited 
(rriB < 14 mag) sample of ~ 7000 galaxies covering galac- 
tic latitudes \b\ > 20° (Marinoni et al. 1999, Giuricin et 
al. 2000). These works presented three different group 
catalogs with groups drawn from the NOG sample us- 
ing a hierarchical grouping algorithm and two different 
variants of a percolation algorithm: one with the linking 
parameters kept constant, and the other with the linking 
parameters scaled with redshift. 

We have cross-referenced the AGC with NOG galax- 
ies with group identifications in the percolation group 
catalog that has redshift dependent linking parameters, 
referred to by Giuricin et al. as the 'P2' group catalog. 
We have combined the P2 group identifications with the 
SCI/SC2 group assignments described in the previous 
section, to generate a single catalog of group assignments. 

Because in most of the sky, the AGC goes signifi- 
cantly deeper than NOG's 14 mag limit, we have also 
developed a modified percolation algorithm that extends 
NOG groups to include nearby AGC galaxies that are 
too dim to be included in the original NOG catalog. The 
AGC is nearly complete down to tub = 15.7 north of 
decl. — —2°, and nearly complete to mg = 15.0 south 
of decl. = —2°, but with a significant number of galax- 
ies dimmer than these completeness levels as well. So 
we have adjusted the linking parameters to represent the 
same level of density contrast as NOG, but with these 
deeper magnitude limits — the relationships between the 
density contrast and the linking parameters are explained 
in Section 2 of Huchra & Geller (1982). We then made 
slight adjustments to the parameters to include the AGC 
galaxies that are even dimmer than the AGC magnitude 
limits. The new linking parameters were then used to 
add AGC galaxies to nearby NOG groups, but not to gen- 
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erate new groups. This procedure is explained in more 
detail by Springob (2006). 

The combined group catalog of SCI/SC2 clusters and 
NOG P2 groups with extensions is used for all SFI++ 
group assignments south of decl. = —2.0°. 

6.2. Voronoi-Delaunay groups 

The catalog described above is useful for numerous 
applications, however it also has some significant draw- 
backs. First, the NOG groups extend to a systemic ve- 
locity of only 6000 km s _1 . The AGC includes a sig- 
nificant number of galaxies beyond this velocity And 
second, the percolation method for generating group as- 
signments has several drawbacks, as pointed out by Mari- 
noni et al. (2002): The technique is insensitive to local 
variations in the density of points, the fine tuning of the 
linking parameters may lead to systematic differences in 
the properties of groups identified at different redshifts 
in the same dataset, and the technique may also identify 
"groups" that are not physically realistic — physically dis- 
tinct concentrations of galaxies connected by long chains 
of galaxies. 

It is precisely such drawbacks inherent in percola- 
tion grouping that motivated Marinoni et al. (2002) 
to propose a more physically motivated group-finding 
algorithm, the Voronoi-Delaunay method (VDM). The 
Voronoi partition, by non-parametrically smoothing 
data, represents an efficient way to measure packing and 
identify as potential group centers the density peaks in 
the galaxy distribution, while the Delaunay mesh, by 
reconstructing the neighborhood relationship between 
galaxies, represents a natural way to assign group mem- 
bers. 

The procedure works as follows: Galaxy group cen- 
ters are identified by peaks in the galaxy density field. 
All galaxies located within a cylindrical volume in red- 
shift space centered on that concentration are assigned 
to the group. According to this strategy, there is no 
need to introduce an arbitrarily chosen global density 
threshold to judge when a given system is formed. In- 
stead, the dimensions of the cylinder are locally scaled 
on the basis of physical considerations, i.e., according 
to the richness- velocity dispersion correlation (e.g., Bah- 
call 1981). The details of the relationship between these 
dimensions and the group richness are derived from semi- 
analytic galaxy formation algorithms applied to N-body 
simulations. The method is described in detail by Mari- 
noni et al. (2002), and was used to identify groups in the 
DEEP2 Galaxy Redshift Survey (Gerke et al. 2005). 

North of decl. = —2.0°, the AGC is nearly complete 
in optical diameter to 1.0' outside the zone of avoid- 
ance (ZOA). We have artificially filled in the ZOA of 
the AGC to eliminate edge effects, and then used C.M.'s 
VDM code to generate VDM group assignments for AGC 
galaxies with a > 1.0' and decl. > —2.0°. 

For the simulation of large scale structure in the ZOA, 
wc first discarded all real galaxies with \b\ < 15°, then 
filled in that region by creating duplicates of all galaxies 
in the region -2.0 < decl. < 40.0°, but shifted by +5 h in 
R.A, and duplicates of all galaxies north of decl. = 40.0°, 
but shifted by 12 h in R.A. Any duplicate galaxies lying 
within \b\ < 15° after having been shifted in R.A. is re- 
tained as a 'synthetic galaxy' in the ZOA. All duplicates 
outside of that region are discarded. 



The VDM algorithm operates on real and synthetic 
galaxies alike, and so several groups are composed ex- 
clusively of synthetic galaxies, or of a combination of 
real and synthetic galaxies. This is necessary for the ini- 
tial group identification. But the final group catalog is 
purged of all synthetic galaxies, and only groups with at 
least two real galaxies have been retained. 

The catalog contains 5423 real galaxies in 1071 groups 
of 2 or more real galaxies. However, only 355 of those 
groups contain 5 or more members. Like the AGC, 
the mean systemic velocities of the vast majority of the 
groups is less than 10,000 km s _1 . However, there are 
some groups with greater distances, including one bi- 
nary pair with a mean systemic velocity of nearly 22,000 
km s _1 . 

6.3. Combined catalog and group statistics 

We emphasize that all template galaxies retain their 
template group assignments regardless of any group as- 
signments in the VDM catalog or previous assignments in 
the SCI/SC2/NOG catalog. For all other SFI++ galax- 
ies, we use SCI/SC2/NOG group assignments for objects 
south of decl. = —2° or systemic velocities of cz < 1000 
km s _1 and VDM group assignments for all objects out- 
side of those limits. 

For this hybrid group catalog, there are 1360 SFI++ 
galaxies in 736 non-template groups. However, the ma- 
jority of these groups include only one SFI++ galaxy. In 
these cases, we still use the group redshift as opposed to 
the galaxy redshift in computing the peculiar velocity, 
in order to remove the effects of any small scale mo- 
tions of the galaxies with respect to the group center of 
mass. Among the 288 groups that include more than 
one SFI++ galaxy only 22 contain more than 5 SFI++ 
galaxies. 

In order to provide the interested reader with an idea 
of the sizes of the groups involved, the data tables pre- 
sented in this work also include the total number of AGC 
galaxies in each group. However, we caution that these 
numbers should not be used for quantitative applica- 
tions. They are only meant as a guide for the reader 
to distinguish between large groups and small groups. 
The numbers of galaxies in the VDM groups does have 
some significance because of the way the VDM grouping 
was done. But because of the inhomogenous nature of 
the AGC, the galaxy counts in the other groups is less 
meaningful. 

Also, as a rough guide to the uncertainty in the groups' 
systemic velocities, we provide a group redshift error 
statistic, derived under the assumption that each galaxy 
redshift serves as an independent measurement of the 
group redshift, neglecting the error introduced by the 
fact that there is some chance that a galaxy's group 
assignment may be erroneous. The error, e cz , on the 
CMB frame velocity, cz group , of a group with Nagc AGC 
galaxies is then given by the sum in quadrature of the 
errors on CMB frame velocities of the individual galaxies: 



Ecz = T lTe (15) 

where ej is the measurement error on the CMB frame ve- 
locity for each individual galaxy in the group. This mea- 
surement error is typically just a few km s _1 for radio 
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rcdshifts, but can be ~ 30 km s _1 for optical redshifts. 
The AGC includes some optical redshifts for which this 
error is unknown. In such cases, we conservatively as- 
sume €i = 40 km s _1 . 

This grouping scheme represents a unique feature of 
SFI++, as past peculiar velocity catalogs have largely 
been restricted to either 'cluster' or 'field' samples, with 
no provision made for groups that have only a handful 
of members with peculiar velocity data. Exploiting such 
grouping data is possible because of the large number of 
objects in SFI++, which exceeds that of almost every 
past peculiar velocity catalog. Willick et al. (1996) made 
the only previous attempt to incorporate galaxy group- 
ing information into the construction of peculiar velocity 
catalogs. 

7. DERIVATION OF PECULIAR VELOCITIES 

7.1. Peculiar velocities of field galaxies 

Peculiar velocity estimates for galaxies not assigned to 
any group are made according to 



Vgal =CZ gal (l-W- Mm ) 



(16) 



which means that the distance to the galaxy (in velocity 
units) can then be expressed as 



T gal — CZg a i ^gal 



(17) 



where cz ga i is the systemic velocity of the galaxy, and dm 
is the difference between the corrected absolute magni- 
tude of the galaxy and the predicted absolute magni- 
tude that one would expect from the template, given 
that galaxy's corrected width measurement. In comput- 
ing dm for individual galaxies, we also make use of the 
same width-dependent morphological correction for ear- 
lier types that is described in Section 5. This treatment 
of the magnitude as a function of the width (rather than 
the other way around) is known as the "forward TF re- 
lation" . 

Peculiar velocity errors are taken as the sum in quadra- 
ture of the absolute magnitude error, velocity width er- 
ror, and an intrinsic TF scatter term. The intrinsic scat- 
ter, in magnitude units, is given by 

e mt = 0.35 - 0.37[log(VF T F) - 2.5] (18) 
as derived in Paper I. 

7.2. Peculiar velocities of groups 

For each group with more than one galaxy in the sam- 
ple, we determine group peculiar velocities by averaging 
the peculiar velocities of each group member, weighted 
by the error on each individual peculiar velocity. That 
is, a group with N$fi++ galaxies with peculiar velocities 
Vi and velocity errors q has velocity 



Nsfi++ 



Vi/Ci 



Nsfi^ 



l/e< 



(19) 



where Vi is computed slightly differently from the method 
explained in Section 7.1. Instead of using the galaxy 
rcdshift as in Equation 16, we use the group redshift 



where cz group is the mean redshift of all galaxies in the 
group, including those AGC galaxies for which we have 
rcdshift data, but not sufficient photometric or spectro- 
scopic data to include in SFI++. 

Using the same mathematical reasoning as in Equation 
15, the group peculiar velocity error is then given by 



SFI+ + 



(21) 



7.3. Malmquist bias 



Vi — cZp roU p(l 10 



0.2dm\ 



(20) 



'Malmquist bias' is the term generally used to refer 
to biases originating from the spatial distribution of ob- 
jects (Malmquist 1924). It arises from the coupling be- 
tween the random distance errors and the density vari- 
ation along the line of sight. Because of these density 
variations, the probability distribution for the distance 
cannot simply be modeled by a gaussian along the line 
of sight, centered on the measured distance. 

There are two types of Malmquist bias that one must 
be concerned with. First, for a given set of selection cri- 
teria, the probability of one's sample including a galaxy 
with a given apparent magnitude will vary with the dis- 
tance to that galaxy, as a result of several factors: 1) 
the luminosity function is not perfectly flat, 2) within a 
given solid angle, there are more galaxies at larger dis- 
tances than smaller distances, and 3) the selection func- 
tion for the galaxies may vary with distance and/or red- 
shift. While this bias is homogeneous across the sky, 
however, the second form of Malmquist bias is inho- 
mogeneous. It arises from the variations in large scale 
structure along the line of sight. Failure to account for 
this type of bias can lead to spurious infall signatures 
onto high density regions. Examples of Malmquist bias 
correction schemes for peculiar velocity samples include 
Freudling et al. 1995 and Park & Park 2006, which cor- 
rect for Malmquist bias in the SFI catalog. 

While it is fairly straightforward to correct for both 
the inhomogeneous Malmquist bias resulting from large 
scale structure (provided one has access to a reconstruc- 
tion of the local density field at hand) and the homoge- 
neous Malmquist bias resulting from the volume effect 
and the luminosity function, one still needs to account 
for the selection criteria of the sample. This is a serious 
problem for our dataset, in that our selection criteria are 
very inhomogeneous. The best we can hope to do is to 
construct ad hoc selection criteria that will mimic the 
observational properties of the catalog, without regard 
to the prior selection criteria that were used to generate 
the sample. 

We have thus adopted a procedure in which we as- 
sume that, whatever the underlying luminosity function 
and sky distribution of sources may be, our selection 
function forces the output catalog to have the particular 
magnitude distribution that is observed at each redshift. 
That is, whatever the underlying distribution of galaxies 
that could be used as TF targets, our selection criteria 
chooses targets according to a probability function that 
must result in the final catalog having precisely the same 
bivariate distribution of redshift and magnitude that our 
catalog has. We can then compute the probability, p(r), 
that a given galaxy is at a distance r by convolving the 
a priori TF-measured distance and associated errorbar 
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with the number of galaxies per magnitude bin at each 
distance (this takes care of the homogeneous Malmquist 
bias) and the large scale density field (this takes care of 
of the inhomogeneous Malmquist bias). 

We emphasize that this approach has been adopted 
by necessity, owing to our inhomogeneous selection cri- 
teria. If the reader wishes to make use of a subsample 
of our data for which a homogeneous selection function 
can be applied, we strongly recommend adopting a dif- 
ferent Malmquist bias correction approach that makes 
use of those selection criteria. Additionally, if one were 
to extract a subset that imposed additional criteria, such 
as excluding galaxies outside of a given distance or red- 
shift range, one would need to impose an additional 
bias correction. So we urge extreme caution in using 
these Malmquist bias-corrected distances. This is in fact 
why we are providing both corrected and uncorrected 
distances — so that the reader will have the ability to use 
his or her own correction procedure if he or she wishes 
to. Having said that, we describe the details of this bias 
correction procedure below. 

7.3.1. Bias correction procedure 

We break up SFI++ into three declination bins, 
decl. < -17.5°, -17.5 < deal. < -2.5°, and decl. > 
—2.5°, I band apparent magnitude bins of width 0.5 mag, 
and redshift bins of width 1000 km s _1 . (The declina- 
tion bins are used because of the variation in complete- 
ness with declination. Within each declination bin, the 
variation in selection criteria from one point on the sky 
to the next is negligible.) We then make the assumption 
that the apparent magnitude distribution within each bin 
of declination and systemic velocity is most likely to be 
the same as the apparent magnitude distribution within 
a corresponding bin of declination and distance. (This is 
likely to be approximately true as, averaged across the 
entire sky, a particular distance bin should be expected 
to have a roughly equal number of galaxies with positive 
and negative peculiar velocities.) 

For each galaxy, we now convolve the probability distri- 
bution pTF{f) of its distance from the TF measurement 
with the probability p m ag(r) of finding a galaxy with its 
apparent magnitude at each possible distance. Further, 
in order to account for the inhomogeneous Malmquist 
bias resulting from large scale structure, we convolve 
these probability distributions with the density field re- 
constructed from the 2MRS (Erdogdu et al. 2006), pro- 
vided to us by Pirin Erdogdu. The probability of the 
corrected distance to the galaxy being rj is then 

p(n) = kiPTF(n)Pmag(n)Plss(n) (22) 

where p; ss (r*j) is the density distribution along the line of 
sight for the galaxy, as given by the 2MRS density field, 
and k\ is a normalization constant, such that E^(rj) = 
1. For a galaxy with TF- measured distance r ga i and 
peculiar velocity error e v , the summation is done over 41 
different distances ranging from r ga i — 2e v to r ga \ + 2e v , 
separated by O.le^. p(ri) is evaluated at each of those 
distances, and the Malmquist bias corrected distance is 
then given by 

rgal-rnalm = ^i=lAlP(n)ri (23) 

with Malmquist bias-corrected peculiar velocity 



Vgal—malm — CZgal T g al—malm 

and distance / peculiar velocity error 



(24) 



^■v—malm 



nJZi=i,4ip(ri)(ri - r gal _ malm ) 2 (25) 

where hi — 1.35 is a coefficient that we include in order 
to correct for the fact that we are evaluating p(r) only 
within the range of r ga i — 2e v to r ga i + 2e v . 

This is the correction that we apply for field galaxies. 
As explained in Paper I, the Malmquist bias for the tem- 
plate clusters is likely to be negligible, so we ignore it. 
However, it should be applied for the field galaxies and 
nontemplate groups. For the groups, we also provide 
Malmquist bias corrected peculiar velocities and errors. 
We use the same relations that are shown in Equations 
19 and 21, but substitute in the Malmquist bias cor- 
rected velocities and errors for the individual galaxies. 
So, a group with N$fi++ galaxies with Malmquist bias 
corrected peculiar velocities f ma /m,i and errors e ma im,i 
would have a group velocity 



^ group — malm 



2~^i—l.NsFl++ ^malm,i / ^malr, 
zJi=l,AT SF/++ l/ e malm,i 



(26) 



where v ma lm,i is computed according to Equation 24, but 
again (as in the case of the uncorrected group velocities) 
using the group redshifts rather than individual galaxy 
redshifts. 

Our Malmquist bias corrected group velocity error also 
follows Equation 21 closely: 



£ group— malm 



y/ N S FI++ 



r++ ' 



(27) 



malm. i 



One final note about this procedure must be made. 
The 2MRS density field only extends out to 20,000 
km s _1 . Additionally, our magnitude and redshift bin- 
ning does not extend beyond 20,000 km s _1 , as the 
number of galaxies becomes too sparse at that point. 
We thus assume a flat magnitude distribution and a flat 
density distribution beyond 20,000 km s _1 . Therefore, 
our Malmquist bias corrections for galaxies with distances 
close to or greater than 20,000 km s _1 tend to represent 
little improvement over the uncorrected distances, and 
should be taken with a large grain of salt. There are, in 
any case, very few SFI++ galaxies with distances in that 
regime. 

8. TF DATA COMPILATION 

Here, we provide the corrected photometric and spec- 
troscopic parameters directly used to compute the pe- 
culiar velocities, as well as the peculiar velocities them- 
selves. We are also making these data available online in 
US National Virtual Observatory (NVO)-compliant ta- 
bles hosted by the Cornell Theory Center 14 . In Table 
2, we present the observational parameters and pecu- 
liar velocities for all galaxies that are not in any of the 
template clusters. This includes galaxies that we iden- 
tify as belonging to nontemplate groups, as described by 

14 see http://arecibo.tc.cornell.edu/hiarchive/sfiplusplus.php 
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Section 6. In Table 3, we present observational param- 
eters and peculiar velocities for all nontcmplatc groups. 
And in Table 4, we present observational parameters for 
all galaxies in the template clusters. Parameters for the 
template clusters themselves can be found in Paper I. 

The format of Table 2 (all galaxies not in the template 
sample) is the same as Table 1 for the first five columns. 
We then add the following parameters: 

Column (6). — Logarithm of the corrected rotational 
velocity width in units of km s _1 , log(WrF)- 

Column (7). — The estimated error on log(WrF), e„. 

Column (8). — The observed I-band magnitude, m b s , 
extrapolated to 8^, before the extinction and face-on 
corrections have been applied that produce the values 
used for TF. (This is the same parameter as Table 1, 
Column 13.) 

Column (9). — The apparent magnitude m;, corrected 
for extinction, and corrected to face-on magnitude. 

Column (10). — The corrected absolute magnitude M/, 
computed assuming that the galaxy is at the distance 
given by its redshift, with H n — 100 km s -1 . 

Column (11). — The estimated error cm on the abso- 
lute magnitude. This is not to be confused with the 
apparent magnitude error given in column (14) of Table 
1, which does not include the error contributions from 
the uncertainty in the inclination, redshift, etc. 

Column (12). — The extinction coefficient 7, in magni- 
tudes, computed according to Equations 2-6. 

Column (13). — The inclination i of the plane of the 
disk to the line of sight, in degrees. 

Column (14). — The velocity of the galaxy in the CMB 
frame, cz ga i, in km s~ x , taken to be the midpoint of 
the spectral line profile, regardless of whether optical or 
radio spectroscopy was used. 

Column (15). — The peculiar velocity v ga i in km s _1 , 
uncorrected for Malmquist bias, as given by Equation 16 
for field galaxies. For galaxies in groups, we provide the 
peculiar velocity estimate before averaging by groups, as 
given by Equation 20. 

Column (16). — The estimated error on the (uncor- 
rected for Malmquist bias) peculiar velocity, e v , in 
km s _1 . 

Column (17). — The (uncorrected for Malmquist bias) 
distance to the galaxy, r ga i, as given by cz ga i — v ga i as 
in Equation 17, in km s _1 , for field galaxies. For galax- 
ies in groups, we provide the distance estimate before 
averaging by groups, cz group - v ga i- 

Column (18). — The Malmquist bias corrected peculiar 
velocity v ga i^ ma i m in km s _1 , as given by Equation 24 
for field galaxies. For galaxies in groups, we provide the 
peculiar velocity estimate before averaging by groups. 

Column (19). — The estimated error on the Malmquist 
bias corrected peculiar velocity, e v - ma i m , in km s _1 . 

Column (20). — The Malmquist bias corrected distance 
to the galaxy, r ga[ 

— malm — CZgal v gal — malmi km S , 

for field galaxies. For galaxies in groups, we provide the 
distance estimate before averaging by groups, cz group — 

^gal—malm • 

Column (21). — Group number that the galaxy is as- 
signed to. For SCI/SC2 clusters, we use the same num- 
bering scheme used for SCI/SC2. For NOG groups, we 
use the NOG P2 numbering scheme, but with 30,000 
added to each group number, so that any group num- 



ber between 30,000 and 40,000 is a NOG group. For 
VDM groups, the numbering scheme begins with group 
number 40,001, so that any group number greater than 
40,000 is a VDM group. If the galaxy does not have a 
group assignment, the number is given as '0'. 

Column (22). — A code indicating whether the spectro- 
scopic data are radio or optical. 'H' for HI spectroscopy 
and 'O' for optical spectroscopy. 

The format of Table 3 (non-template groups) is as fol- 
lows: 

Column (1). — Group number, using the same number- 
ing scheme as Table 2, column (21). 

Columns (2) and (3). — Right ascension (in hh mm ss.s 
format) and declination (in dd mm ss format) in J2000.0 
epoch of the group center. 

Column (4). — The CMB frame velocity, cz group , in 
km s _1 , taken to be the average of the CMB velocities 
of the constituent galaxies, regardless of their inclusion 
in SFI++. 

Column (5). — The 'error' on the CMB frame velocity, 
e cz , in km s _1 , computed according to Equation 15. 

Column (6). — The number of galaxies in the group, 
Nsfi++, for which we have data in SFI++. 

Column (7). — The number of AGC galaxies in the 
group, N AGC . 

Column (8). — The group peculiar velocity v group , in 
km s _1 , uncorrected for Malmquist bias, computed ac- 
cording to Equation 19. 

Column (9). — The estimated error on the (uncorrected 
for Malmquist bias) group peculiar velocity, e groupi in 
km s _1 , computed according to Equation 21. 

Column (10). — The (uncorrected for Malmquist bias) 
distance to the group, r group , as given by cz group - v gr0 u P , 
in km s _1 . 

Column (11). — The Malmquist bias corrected group 
peculiar velocity v group -. ma i mi in km s _1 , uncorrected 
for Malmquist bias, computed according to Equation 26. 

Column (12). — The estimated error on the Malmquist 
bias corrected group peculiar velocity, e group - ma i m , in 
km s _1 , computed according to Equation 27. 

Column (13). — The Malmquist bias corrected distance 
to the group, r group — n i a i m cz group u group — nia i mj in 
km s _1 . 

The format of Table 4 (galaxies in the template sam- 
ple) follows that of Table 2, except that we omit columns 
15-22 and add columns (23), the angular separation, 8, 
between the galaxy and the cluster center, in units of ar- 
cminutes, and (24), the template cluster name, using the 
same naming scheme as in Paper I. We have also added 
column (25), which contains both the spectroscopic data 
code from Table 2, column (22), and a second code which 
characterizes membership status in the cluster. As in 
G97, code 'c' signifies a bona fide cluster member (the 
in sample), while code 'g' indicates that a firm mem- 
bership assignment cannot be made (the in+ sample). 
One other difference between Tables 2 and 4 is that, fol- 
lowing G97, we use the cluster redshift to compute the 
absolute magnitude of in galaxies in column (10). For 
in+ galaxies, we use the galaxy redshift, as in Table 2. 

9. DATASET CHARACTERISTICS 

The sky distribution of all SFI++ galaxies is shown 
in Figure 2. As with the S05 HI archive, very few ob- 
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Fig. 2. — Aitoff equal area projection of the sky distribution of 
SFI++ galaxies. The plot is centered at R.A. = 6 h . The thick lines 
trace the galactic latitudes b = -20°, b = 0°, and b = +20°. The 
figure is color-coded by CMB frame redshift, such that cz ga i < 1000 
km s — 1 galaxies are blue, 1000 < cz ga i < 3000 km s — 1 galaxies 
are cyan, 3000 < cz ga i < 5000 km s — 1 galaxies are green, 5000 < 
cz gal < 7000 km s _1 galaxies are yellow, 7000 < cz ga i < 9000 
km s — 1 galaxies are red, and cz ga i > 9000 km s — 1 galaxies are 
magenta. 

jects can be found close to the Galactic plane. SFI++ 
does, however, include a far greater share of galaxies 
in the southern hemisphere than the HI archive; how- 
ever there is still a deficiency of galaxies in the range 
— 17.5 < decl. < —2.5°, owing to the fact that all of our 
observing programs targeted known galaxies in published 
catalogs, and there is no catalog of comparable depth to 
the UGC or ESO in this declination range. We also note 
that, due to the large amount of HI data we have in 
the Arecibo declination range, there is better coverage 
there (—2 < decl. < +38°) than in any other part of the 
sky. In Figure 3, we show the sky distribution of the tem- 
plate clusters. The clusters are broadly spread across the 
sky, though there is a greater concentration in both the 
Pisces-Perseus Supercluster and the supergalactic plane. 

Figure 4 shows the distributions of apparent and ab- 
solute magnitudes, widths, radial velocities, and TF dis- 
tances, both corrected and uncorrected for Malmquist 
bias. Because of the redshift-dependent diameter limit 
of the SFI, specifically designed to equalize the number 
of objects across a large range of distances, the radial 
velocity distributions of Figure 4d is somewhat "flatter" 
than the corresponding radial velocity distribution of the 
HI archive shown in S05 Figure 4. As in that figure, the 
redshift distribution peaks at ~ 5000 km s _1 , which 
corresponds to the mean recessional velocity of the dens- 
est concentration of the Pisces- Perseus Supercluster. Be- 
cause the TF distances for individual galaxies have errors 
of ~ 15%, that redshift peak is even more smoothed out 
in the TF distance histogram, though moreso for the his- 
togram of uncorrected distances (Fig. 4e) than the his- 
togram of Malmquist bias corrected distances (Fig. 4f). 

The distributions of the peculiar velocities of individ- 
ual galaxies, nontemplate groups, and template clusters 
are shown in Figure 5. (In the case of individual galaxies 
and nontemplate groups, we provide the distributions of 



Fig. 3. — Same as Figure 2, but for SFH — f- template clusters. 
We use the same coloring scheme as in Figure 2, color-coded by 
the mean cluster redshift. 
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Fig. 4. — Distribution of (a) Observed I band magnitude, in 
bins of width 0.25 mag, (b) I band absolute magnitude, in bins of 
width 0.25 mag, (c) logarithm of velocity widths, in bins of width 
0.025 dex, (d) recessional velocity in CMB frame, in bins of width 
500 km s _1 , (e) TF distance (uncorrected for Malmquist bias), 
in bins of width 500 km s _1 , and (f) Malmquist bias corrected 
TF distance, in bins of width 500km s _1 , for all galaxies in the 
SFI++. 

peculiar velocities in both the Malmquist bias corrected 
and uncorrected cases.) We emphasize that very few 
galaxies have real peculiar velocities of magnitude greater 
than 1000 km s _1 . In fact, few galaxies outside of clus- 
ters have real peculiar velocities of magnitude greater 
than 500 km s _1 . Thus, the broad distribution of pe- 
culiar velocities for individual galaxies and small groups 
is a consequence of the scatter in the TF relation — some 
galaxies are intrinsically offset from the TF relation, and 
some simply have unusually large magnitude or width 
errors. 

The peculiar velocities of template clusters fall in a 
much narrower range, both because of the reduced dis- 
tance errors associated with these clusters, and the fact 
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Fig. 5. — Top Peculiar velocity distribution of all nontcm- 
platc groups (hatched) and all field galaxies (unhatched) before the 
Malmquist bias correction; middle peculiar velocity distribution of 
all nontcmplatc groups (hatched) and all field galaxies (unhatched) 
after the Malmquist bias correction; bottom peculiar velocity dis- 
tribution of the template clusters. In all three plots, the bins width 
is 100 km s" 1 . 

that larger groups and clusters tend to have smaller pe- 
culiar velocities than do field galaxies or small groups. 
The dispersion of template cluster peculiar velocities is 
discussed in Paper I. 

The TF relation itself is plotted in Figure 6, superim- 
posed against Paper I's template relation. All nontcm- 
plate galaxies arc included in this figure, with the mag- 
nitudes for earlier type galaxies corrected by the width- 
dependent morphological corrections of Section 5. Be- 
cause this plot does not account for the Malmquist bias 
correction, selection effects should be apparent, and we 
expect there to be some deviation from the template re- 
lation. However, there are some extreme outliers — many 
of which have unusually large width errors. Most, but 
not all, of these galaxies with large width errors have 
small inclination angles, so the corrections to edge-on 
widths are large. Some simply have larger-than-average 
measurement errors — these are far more likely to be op- 
tical widths than radio widths. There do remain some 
outliers, however, that do not have unusually large mag- 
nitude or width errors. It is possible that in such cases, 
we have underestimated the uncertainties. One compli- 
cating factor is that, since the TF relation depends on 
morphological type, ambiguity in morphological classifi- 
cation will inevitably introduce its own sources of error, 
though this error is extremely difficult to quantify. This 
problem is more significant for this dataset than in most 
of our group's previous TF studies, because SFI-I — h in- 
cludes more galaxies of smaller optical extent, which are 
more difficult to classify. In any case, the extreme out- 
liers from the TF relation would likely be excluded for 
most applications of the catalog, but we leave it up to 
the reader to decide on an appropriate cutoff. 

10. SUMMARY 

We have constructed a TF catalog, the SFI++, includ- 
ing 4861 field and cluster galaxies, primarily based on 
observations made by our group and collaborators over 
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Fig. 6. — Top TF relation for all nontcmplatc galaxies. Bottom 
Same points plotted, but with magnitude and width crrorbars in- 
cluded. The solid line represents the template relation from Paper 
I, which wc provide in Equation 14. The width-dependent morpho- 
logical corrections of Section 5 have been applied to the magnitudes 
plotted here, unlike the magnitudes in Table 2 and 4. 

the course of the last ~ 25 years, but also including some 
data from the literature that has been processed in an 
identical manner. While most of the spectroscopic data 
have already been published, all of the data have now 
been reprocessed using new correction schemes. And a 
significant fraction of the photometric data are provided 
here for the first time. We have reexamined the choice 
of clusters used in the construction of a TF template, 
as well as the group assignments of all galaxies in the 
vicinities of those clusters. A new template TF relation 
has been derived, as presented in Paper I. 

We have also used the 2MRS density field to correct our 
measured distances and peculiar velocities for Malmquist 
bias. The corrections require us to construct ad hoc se- 
lection criteria for our sample, and we strongly advise 
that if one is to make use of a subset of our data that 
involves imposing additional selection constraints, then 
one should not use the same bias corrections presented 
here. Additionally, our Malmquist bias corrections are 
likely to be of limited value for galaxies with distances 
close to or greater than 20,000 km s _1 . (Though very 
few galaxies in our sample are found in this regime any- 
way.) 

The complete TF catalog is by far the largest such 
peculiar velocity catalog with ~ 15% distance errors to 
date. SFI++ is also unique among peculiar velocity cat- 
alogs in that it includes a mix of cluster and field galax- 
ies, as well as group identifications for galaxies in loose 
groups. 

The detailed spectroscopic datasets, including parame- 
ters not directly used in the computation of TF distances, 
have been presented in S05 and Catinclla, Haynes, & Gio- 
vanclli (2005), but we have provided the detailed photo- 
metric parameters for most of the galaxies in our sam- 
ple here, including galaxies for which we do not have 
high quality width measurements, which are therefore 
excluded from the TF sample. The distribution of pa- 
rameters is very similar to those provided in subsamplcs 
of the data published by our group before we made re- 
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cent changes to the data processing algorithms. The en- 
tirety of the TF dataset is available online through NVO- 
compliant tables hosted by the Cornell Theory Center. 
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